Probing quantum-mechanical level repulsion in disordered systems by means of 
time-resolved selectively-excited resonance fluorescence 
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We argue that the time-resolved spectrum of selectively-excited resonance fluorescence at low 
temperature provides a tool for probing the quantum-mechanical level repulsion in the Lifshits tail 
of the electronic density of states in a wide variety of disordered materials. The technique, based 
on detecting the fast growth of a fluorescence peak that is red-shifted relative to the excitation 
frequency, is demonstrated explicitly by simulations on linear Frenkel exciton chains. 
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Anderson localization is a key concept in understand- 
ing the transport properties and optical dynamics of 
quasiparticles in disordered systems in general [H, 0, E| 
and in a large variety of low-dimensional materials in 
particular. Examples of current interest are conjugated 
polymers [1| , molecular J-aggregates [!, 0] , semiconduc- 
tor quantum- wells 0, H| and quantum wires 0] , as well 
as biological antenna complexes [lfj, . Localization re- 
sults in the appearance of a Lifshits tail in the density of 
states (DOS) below the bare quasiparticle band (lil. fl3j|. 
At low temperature, the states in the tail determine the 
system's transport and optical properties. Of particular 
importance for the dynamics and the physical properties 
is the spatial overlap between these states. Two situ- 
ations can be distinguished: states that do not overlap 
can be infmitesimally close in energy, while states that 
do overlap undergo quantum-mechanical level repulsion. 
Interestingly, this repulsion does not manifest itself in 
the overall level statistics, which, due to the localization, 
is of Poisson nature. The local Wigner-Dyson statistics, 
caused by the repulsion, turns out to be hidden under 



the global distribution of energies (14|. Still, due to the 
characteristic energy scale associated with the level re- 
pulsion, this phenomenon does affect global properties, 
such as transport [lij]. Moreover, observing changes in 
the level statistics allows one to detect the localization- 
delocalization transition or mobility edges 0]. 

Thus, it is of general interest to have experimen- 
tal probes for the level statistics in disordered sys- 
tems. Recently, it has been shown that near-field spec- 



troscopy flil. 171 and time-resolved resonant Rayleigh 
scattering [l8j may be used to uncover the statistics 
of localized Wannier excitons in disordered quantum 
wells [l(| and wires [T3|- In this Letter, we argue that 
low-temperature time-resolved selectively excited fluores- 
cence from the Lifshits tail provides an alternative tool 
for probing the level repulsion. This method is based 
on the fact that downward relaxation between spatially 
overlapping states dominates the early-time rise of a fluo- 
rescence peak that is red-shifted relative to the excitation 



frequency. The red shift is related to the level repul- 
sion. We will demonstrate this by simulations on a one- 
dimensional (ID) Frenkel exciton model with on-site dis- 
order. The key ingredients of the model - level repulsion 
and scattering rates proportional to the phonon spectral 
density times the overlap integral of the site probabilities 
of the two exciton states involved [6|_- are also charac- 
teristic for Wannier exciton systems fH IH lil. Il9j]. Hence, 
the method has a wide range of applicability. 

Our model consists of an open chain of N optically 
active two-level units (monomers) with parallel transi- 
tion dipoles, coupled to each other by dipole-dipole in- 
teractions. The chain's optical excitations are Frenkel 
excitons described by a Hamiltonian matrix whose di- 
agonal elements are the monomer excitation energies, 
H nn = £„ , which are random uncorrelated Gaussian vari- 
ables with zero mean and standard deviation A (i.e., 



(e n ) = and (e„e TO ) = A 2 8 n 



ith 



denoting the 



average over the disorder realizations). The off-diagonal 
matrix elements are non-random dipole-dipole interac- 
tions: H nm = —J/\n — m\ 3 (J nn = 0), with J > 
representing the magnitude of the nearest-neighbor cou- 
pling. The exciton cigenencrgies E v and eigenfunctions 
(f vn follow from diagonalizing H nm . 

The states in the Lifshits tail, which dominate the op- 
tical absorption of the disordered chain, are localized on 
segments of typical size N* . Some of these states can be 
joined in local exciton manifolds of a few states that are 
similar to the eigenstates of an ideal chain of size equal to 
the segment length [lilliol]. The states in a local manifold 
undergo level repulsion, while between different segments 
energy separations may be arbitrarily small. Because the 
variation of exciton energies over different segments is 
larger than the level repulsion within segments [3], the 
latter remains hidden in techniques that probe global av- 
erages like, e.g., the absorption spectrum. 

The typical situation is illustrated in Fig. [TJ The two 
states to the right form a local manifold, with the lower 
state having no node and collecting most of the oscilla- 
tor strength of the monomers within the segment. We 
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FIG. 1: Cartoon of the local DOS tail states, illustrating the 
possible channels of the exciton relaxation at zero temper- 
ature after excitation within a narrow window. The solid 
(dashed) arrow shows the intra (inter) segment relaxation 
with the rate Wsv (Wjg's); while the wavy arrows denote 
the radiative decay with the typical rate 7g = 70-/V* and 
7p ~ O.I70-/V* (N* is the typical segment length). 



will refer to this as an S type state. The higher state 
has one node (P type) and has a much smaller (but 
nonzero) oscillator strength; thus, it can be photo-excited 
but its radiative decay is slow compared to that of the 
S state. The probability overlap between the two states, 



measured by the integral = ^Z n =i ^\nAni is 011 tnc 
order of 1/N seg , with N scg ~ N* the segment (localiza- 
tion) size. The other two states in Fig. [T] are the lowest 
(S type) states of other localization segments. For linear 
systems, the overlap between states of adjacent segments 
typically is two orders of magnitude smaller than the 
overlap within a manifold 14j; for distant segments the 
overlap is even much smaller. These overlap properties 
play a crucial role in the intraband exciton relaxation. 

To describe this relaxation, we use the Pauli master 
equation for the populations P v of the exciton states: 



P v = R V F V - 7 „P„ 



N 

£ 

^=1 



(1) 



where i?„ is a source term, specified below, F v = 
(52n=i Vvn) 2 is the dimensionless oscillator strength of 
the vth. exciton state, 7„ = 70-Fl/ is its spontaneous 
emission rate (70 being the spontaneous emission con- 
stant of a monomer), and Wa V is the scattering rate 
from the exciton state \v) to the state \fx), which re- 
sults from weak exciton-vibration coupling. The lat- 
ter may be obtained through Fermi's golden rule Q: 
W^u = Ifj, v SQu^vl) Gfanv), where 1^ is the overlap 
integral defined above and SQui^]) is the phonon spec- 
tral density at the energy u)^ v — — E v . Furthermore, 
G(u>) = n(u) if uj > and G(u>) = 1 + n(-w) if u < 0, 
withn(w) = [exp(w/T)-l] _1 , the mean thermal occupa- 
tion number of a phonon mode of energy u> (H = ks = 1)- 
We consider an experiment in which a short narrow- 
band laser pulse is used to selectively excite states in 
a small part of the absorption band and set R v (t) = 
RoS(t)S(Ei — E v ) with Ei being the laser frequency. We 



are interested in the time-resolved fluorescence spectrum 
at low temperature (T — 0). The spectrum is given by 



I(E,t) = ±{j2l»P4t) S(E-E U ] 



N 



N\ 



(2) 



At T — 0, only states with E v < Ei are relevant. Quite 
generally, the spectrum consists of a narrow resonant 
peak, resulting from fluorescence from initially excited 
states and a red-shifted broader feature, arising from 
states that are populated by downward relaxation from 
the initially excited ones. 

If the laser frequency lies in the main part of the ab- 
sorption band, not too far in the blue wing, the initially 
excited states are either of the S or of P type. Excita- 
tion of higher excited states, which extend over various 
segments, may be neglected. Both S and P states con- 
tribute to the narrow resonant fluorescence, with most of 
the intensity coming from the S type states, because they 
carry more oscillator strength. On the other hand, the 
P type states have a much higher probability to relax to 
their S partner before they decay radiatively. The S type 
states can only relax to states of other segments, which is 
a slow process as compared to the intra-manifold relax- 
ation (Wsp > Ws's)- Therefore, at times t < 1/Wsp, 
any red-shifted fluorescence comes from S type states 
that have been excited via intra-segment relaxation from 
P type ones. Thus, the shape and position of the red- 
shifted fluorescence band reflects the level repulsion, i.e., 
its intensity tends to zero when approaching the laser 
frequency. At longer times, this repulsive gap may be 
partially filled due to inter-segment hops. 

We now analyze in more detail the early-time shape of 
the red-shifted band I Ied (E,t). For t < 1/7^ + 1/Wsp, 
one iteration of the master equation gives 

I Ied (E,t) = ^{^^W vll F IJ ,5{E ll -E l )5{E-E v 

fl.U 



= f^SiEt-E^^V, 
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x 5{E )1 — Ei) S(E — E v ) 



(3) 



where the sum is restricted to states which occur in SP 
doublets, as described above. Thus, /i (y) in the sum 
label P (S) type states. Because Isp ~ l/JV seg , while 
the lower state in each segment has a superradiant decay 
rate 7s « 7o^Vscg, we have 7„ I Vfl = 070, with a being 
a constant of the order unity. If we now replace F^ by 
the average oscillator strength Fp(Ei) for P type states 
at the energy Ei, we arrive at 

I ied (E,t) = Fp ( El ) S(Et - E) P SF (E h E) . (4) 

Here, P SP (E h E) = { £^ 8(E„ - E t ) 5{E - E v )) is the 
probability distribution for the energy of an S state under 
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the condition that the same localization segment contains 
a P state at the energy E\. This result shows that the 
early-time lineshape is the product of the phonon spec- 
tral density and the conditional spacing distribution, thus 
containing detailed information about the repulsive level 
statistics. 

To confirm the above analytical result, we have per- 
formed numerical simulations of the spectrum Eq. (|2|) for 
ID Frenkel chains of length N = 1000 at T = 0. In all 
calculations we considered system parameters that are 
appropriate for linear aggregates of the dye pseudoiso- 
cyanine (PIC) with counter ion F~: 70 = 2.7 x 10 8 s _1 = 
1.5 x 10~ 5 J, J = 600 cm" 1 , and A = 0.25 J, giving a 
typical localization size N* — 28. For the spectral den- 
sity we used the Debye model: S(u>) = Wq(lu/J) 3 , with 
Wo = 15. 9 J. This set of parameters has been used suc- 
cessfully to analyze and fit experiments on aggregates of 
PIC-F 0. As at T — only states below the excita- 
tion energy are relevant, we used the Lanczos method to 
calculate the subset of only these eigenstates. The fluo- 
rescence spectra were obtained after averaging over 10 5 
realizations of the on-site energies. 



was chosen such that it pumps primarily P type states. 
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FIG. 2: Calculated energy distributions of S type states 
(solid curve), P type states (dash-dotted curve), absorption 
spectrum (dashed curve), and the conditional probability 
P SP (£,,£) (filled circles). 

First, we analyzed numerically various statistical prop- 
erties of the exciton levels. The results are summarized 
in Fig. [5] Shown are the distribution of energies of the 
S and P type states, Ps(E) and Pp(E), respectively, the 
absorption spectrum, J a b s (-E), and the conditional prob- 
ability P SP (Ei,E), with Ei = -2.41J. The selection of 
doublets of S and P type states was done according to the 
method described in Ref. 14 2l[ . Clearly, the low-energy 
half of the absorption peak is dominated by S type states, 
while the high-energy wing contains also important con- 
tributions from P type states. Both distributions overlap, 
which is the reason that the local energy structure is not 
seen in the absorption peak. The conditional probability 
clearly demonstrates the level repulsion, as it drops to 
zero at E = Ei. In our simulations the laser energy Ei 
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FIG. 3: Zero-temperature time-resolved fluorescence spectra 
at different instants. The intensities of the red-shifted bands 
are normalized to the same peak value. The inset shows the 
unnormalized spectra. 

Figure [3] presents our numerical results for the time- 
resolved fluorescence spectrum after selective excitation 
within a window [—2. 415 J, —2. 405 J]. Indeed, the spec- 
trum reveals a sharp peak at the excitation energy and a 
clearly separated red-shifted (R) band, indicative of level 
repulsion. The kinetics of this band can be understood 
from estimates of the radiative decay and relaxation 
rates. Typically, 7^ = -y Q N* w 107 P w 4 x 10~ 4 J, while 
W SP « 10 2 IT S 'S ~ W {A*/J) 3 I SP re 6 x 10~ 5 J. Here 
we used the typical energy separation A* = A/V N* and 
J BP a 1/N* (N* = 28). Thus, W S 's < 7p < W S p < 7s . 
implying that inter-segment hops are irrelevant on the 
radiative time scale. As a consequence, the observed 
blue-shift of the R band upon increasing time, cannot 
be the result of filling the repulsive gap by such hops. 
Instead, this kinetics results from the fact that the oscil- 
lator strength per state peaks at lower energies than the 
absorption maximum [6j. Thus, the higher-energy states 
in the R band decay slower, causing the blue-shift. 

Finally, Fig. [¥] displays the numerically obtained early- 
time fluorescence spectra, for t = O.OlWgjt, 0.02WgJ T, 
and 0.03W^> . The spectra are normalized by dividing by 
the time t at which they were taken. Clearly, the result- 
ing curves cannot be distinguished. Moreover, they agree 
almost perfectly with the function (Ei — E) 3 Psp (Ei , E) , 
plotted as open circles. This confirms the analytical ex- 
pression (j4|). Without showing explicit results, we note 
that the same quality of agreement has been found for 
other model parameters and spectral densities. 

In summary, we have shown that time-resolved selec- 
tive excitation of fluorescence in the Lifshits tail of low- 
temperature disordered systems, yields a probe for the 
quantum level repulsion and local level statistics. Shortly 
after the laser pulse, the shape of the red-shifted fluo- 
rescence band arising in this technique is the product 
of the level spacing distribution within localization seg- 
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FIG. 4: Simulated fluorescence spectra at three early in- 
stances of time, normalized by dividing by t and compared 
to the function (E t — E) 3 Psp(Ei, E) (open circles). The inset 
shows the unnormalized spectra. 

merits and the phonon spectral density. To extract the 
spacing distribution, the spectral density must be mea- 
sured. The widely adapted scheme for this also employs 
low-temperature selectively excited fluorescence, but un- 
der steady-state conditions. It is assumed that the red- 
shifted band directly reflects the phonon spectral den- 
sity [22| . However, for multichromophore systems with 
a fast intra-band relaxation, this approach is not ap- 
propriate, because a relaxation-induced red-shifted fea- 
ture, which also arises under steady-state conditions [23| . 
is superimposed with the phonon side-band. The two 
contributions may be separated by measuring the time- 
dependent spectrum. Thus, time resolution is an essen- 
tial aspect of the scheme proposed. Of course, the rela- 
tion between pulse width (selective excitation) and dura- 
tion (time-resolution) limits application of our approach 
to systems for which the absorption band is wide com- 
pared to the intraband relaxation rate. For the example 
considered here, this criterion is easily met. 

To conclude we stress that our method has a wider ap- 
plicability than the ID Frenkel exciton system we mod- 
eled. The reason is that the key ingredients - the nature 
of the states in the Lifshits tail and the importance of 
spatial overlap for their phonon-induced relaxation - are 
shared by a large variety of systems d, 0, H, H> Ell • 
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